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Abstract. We suggest that a magnetic-field-induced Peierls instability accounts 
for the recent experiment of Zhang el al. in which unexpected quantum Hall 
plateaus were observed at high magnetic fields in graphene on a substrate. This 
Peierls instability leads to an out-of-plane lattice distortion resulting in a charge 
density wave (CDW) on sublattices A and B of the graphene honeycomb lattice. 
We also discuss alternative microscopic scenarios proposed in the literature and 
leading to a similar CDW ground state in graphene. 

1 Introduction 

The recent interest in graphene was triggered by the observation of an unusual "relativistic" 
quantum Hall (QH) effect by two groups |l|2j . This effect is a direct consequence of the peculiar 
density of states of graphene close to half-filling [3]. Shortly after these pioneering experiments, 
Zhang et al. [I] explored the regime of even larger magnetic helds, up to 45 T, and observed 
extra plateaus in the Hall conductivity a xy . In contrast to the relativistic quantum Hall effect 
observed at lower magnetic fields, these extra plateaus cannot be explained by considering 
non-interacting electrons. In the present paper, we review the recent proposal of a magnetic- 
field-induced Peierls instability for graphene on a substrate [5], which gives an explanation to 
the extra quantum Hall plateaus observed at high magnetic field. The Peierls instability leads 
to an out-of-plane lattice distortion resulting in a charge density wave (CDW) on sublattices A 
and B of the graphene honeycomb lattice. Other mechanisms leading to such a CDW ground 
state in graphene have been proposed in the literature and are discussed in the present paper 
|6l7l8l9j . The originality of the magnetic- field- induced Peierls instability scenario 5j is to rely 
on electron-phonon coupling via the substrate on which the graphene flake lies, rather than on 
Coulomb interactions between electrons 6 7 8 9. 

2 Quantum Hall effect in graphene: non-interacting electrons 

The graphene crystal is a honeycomb lattice of carbon atoms: a two-dimensional triangular 
Bravais lattice with a basis of two atoms, usually referred to as A and B. The distance between 
nearest neighbor atoms is a — 0.14 nm and the lattice constant is a\/3. Experimentally, graphene 
sheets of area A ~ (3—10 /im) 2 are deposited on amorphous Si0 2 /Si substrate. Applying a gate 
voltage V g via the substrate allows one to control the electronic filling of the graphene bands, 
which amounts to doping. The number of induced electronic charges is given by N c = V g C g /e 
where the capacitance per unit area can be estimated as C g /A ~ e/d 1.2 x 10~ 4 F/m 2 , where 
— e < is the electron charge, e = eoe r (with e r m 4) is the silicon oxide dielectric constant and 
its thickness is d ~ 300 nm [1 12] . 
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The electronic properties of graphene are described by a standard nearest-neighbor tight- 
binding model [10] with hopping amplitude t ss 3 eV [IT] between 2p z carbon orbitals. Each 
carbon atom contributes one conduction electron. There are 2N p electrons in the sample under 
zero gate voltage, where N p is the number of unit cells in the sample. The first Brillouin zone 
is hexagonal. Among its six corners, only two are inequivalent and called K and K' . We choose 
K = 47r/(3v / 3a)u K and K' = — K. The resulting band structure features the merging of the 
conduction and valence band at precisely these two points: graphene is a two-valley (K and 
K') zero-gap semiconductor. Near these so-called Dirac points, the electrons behave as charged 
massless chiral fermions and the valence and conduction bands have the shape of a diabolo. 
The relativistic-like dispersion relation of electrons is £k = ±Hvf\\s\, with a Fermi velocity 
vp = 3at/2h w 10 6 m/s about 300 times smaller than the velocity of light c. This results in 
a linear density of states close to the Dirac point. In undoped graphene (V g — 0), the valence 
band is filled and the conduction band is empty: the chemical potential is right at the Dirac 
points. Changing the gate voltage from zero allows one to fill the conduction band or empty 
the valence band. 

In order to discuss the quantum Hall effect in graphene, we briefly review the Landau 
problem of a graphene electron in a perpendicular magnetic field [12] . The magnetic field is 
assumed weak enough to neglect any magnetic flux commensuration effect (for a study of the 
Hofstadter butterfly in graphene, see |13j ) and we use the continuum approximation in which the 
tight-binding model reduces to the two-dimensional massless Dirac Hamiltonian. The Landau 
levels (LL) are 

e n = sgn(n) v /jn|fiw c , (1) 



where the "cyclotron" frequency is u> c = vpy / '2eB±/h, and the LL index n is an integer. The 
unusual LL structure with a square root dependence in the perpendicular magnetic field is a 
direct consequence of the linear density of states close to the Dirac point. Each LL is AN^ times 
degenerate, where the number of flux quanta across the sample = eB±A/ \2ttK) gives the 
orbital degeneracy and the factor 4 accounts for twofold spin 1/2 and twofold valley degeneracy. 

Let v = Nc/Nff, = CgVg/eNff, be the filling factor. When the gate voltage is zero, so is the 
filling factor and the central Landau level (CLL, n = 0) is half- filled as a result of particle- 
hole symmetry. In addition, each LL is AN^ times degenerate. These two facts alone show 
that a plateau at ve 2 /h occurs in the Hall conductivity a xy each time the filling factor is 
close to sgn(n)(4|n| + 2) = ±2; ±6; ±10; .... In order to have extra quantum Hall plateaus, 
one has to provide mechanisms for lifting the spin and valley degeneracy. In the absence of 
interactions between electrons, the only such mechanism is the Zeeman effect, which lifts the 
spin degeneracy, but does not affect the twofold valley degeneracy. The Zeeman effect, which 
depends on the total magnetic field Btot, splits each LL and leads to a gap Az = g^s-Btot, 
where [Ib — eh/2m is the Bohr magneton, m is the bare electron mass, and g* « 2 is the 
effective ^-factor [3] . If the Zeeman gap is larger than the LL width, one expects extra quantum 
Hall plateaus at v = An — 0; ±4; ±8; .... Therefore in the absence of interactions between 
electrons, the only expected quantum Hall plateaus should occur at every even filling factor 
v = 2n. In the experiment of Zhang et ai, plateaus at v = ±1 were observed at large magnetic 
field, suggesting that interactions between electrons play a role [4]. However, none of the other 
odd plateaus v = ±3; ±5; . . . were observed. It is very simple to understand these facts if one 
assumes that the electrons have different on-site energies on the two carbon atoms in the unit 
cell, just as for boron nitride, which also has a honeycomb structure |14J . In that case the Dirac 
electrons become massive Q3] and the LL are almost unchanged except for the CLL which is 
not valley degenerate anymore [15]: this is enough to understand the filling factors at which 
QH plateaus were observed. It now remains to find a microscopic mechanism leading to the 
generation of this mass. In the following section, we show that in the presence of a magnetic field 
phonon-mediated (indirect) interactions between electrons lead to a CDW ground state. The 
latter features different electronic charge densities on the two sublattices providing different 
on-site energies for A and B carbon atoms, and therefore leading to a mass for the Dirac 
electrons. 
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3 Magnetic-field-induced Peierls instability 

Undoped graphene is a very peculiar two-dimensional metal: it has a Fermi surface consisting of 
only two points (K, K'), which means that its density of states at the Fermi energy is vanishingly 
small. Such a Fermi surface features perfect nesting at wavevector Q = and K, as testified by 
the Kohn anomalies found in the spectrum of its optical in-plane phonons [16] . This means that 
graphene is a quite unstable metal with respect to charge density modulations at wavevector 
Q and therefore, although two-dimensional, a good candidate for a Peierls instability [TJj. 
Electrons in free-standing graphene are strongly coupled to the in-plane optical phonons but 
are not coupled to the out-of-plane ZO phonon [16j . Such a coupling is also absent in graphite 
because each graphene layer is in a symmetric environment with respect to its neighboring 
layers. However for graphene on a substrate, this is not the case: as the substrate (amorphous 
SiC>2) is different from the superstrate (air), the coupling of electrons to the ZO phonons is 
not forbidden by symmetry. Such an electron-phonon coupling is actually provided by the 
interactions of the graphene sheet with the substrate, as we show below. The ZO phonons 
correspond to an out-of-plane vibration where the A and B atoms in the unit cell move out of 
phase. If this vibration is frozen-in, every second carbon atom is closer to the substrate than 
its partner in the unit cell: the honeycomb lattice is no longer flat but crinkled as an egg box. 
In that case, the A and B carbon atoms do not have the same on-site energy anymore, which 
means that the Dirac electrons have acquired a finite mass. In the following, we quantitatively 
develop this idea. 

We now consider a spontaneous out-of-plane lattice distortion which - in the presence of the 
substrate - breaks the inversion symmetry of the honeycomb lattice and provides a mechanism 
for lifting the valley degeneracy in the CLL. Assume that the A (resp. B) sublattice moves away 
(resp. towards) the substrate by a distance r/. Electrons are still described by a honeycomb 
nearest-neighbor tight-binding model, however the two atoms in the basis now have different 
on-site energies. The energy on atom A/B is called ±M. Close to the Dirac points and including 
the Zeeman effect, the LL read [15] 



where a — ±1 is the valley index corresponding to the Dirac points aK, and a = ±1 is the spin 
projection along the magnetic field direction. In terms of the low-energy effective theory, the 
distortion means that the Dirac electrons spontaneously acquire a finite mass M/v 2 F . The on- 
site energy difference lifts the valley degeneracy for the CLL only, with a valley gap A v = 2M. 
Note that the effect of a nonzero on-site energy M on each n ^ LL is very weak, of order 
M 2 /HvpeB± <C 1 for a typical magnetic field. The reason is that in every LL, apart from the 
CLL, each quantum state - independently of its valley index - has an equal weight on the two 
sublattices. This is, however, not the case in the CLL: electrons at the K (resp. K') point only 
reside on the A (resp. B) sublattice. This explains why breaking the inversion symmetry of 
the honeycomb lattice (by introducing different on-site energies on the two sublattices) lifts the 
valley degeneracy of the CLL only. 

Such a lattice distortion spontaneously occurs because it lowers the total energy, just as in 
Peierls's mechanism [17] except for the magnetic field playing an essential role here and for the 
crystal being two rather than one dimensional. Assume that the last partially filled LL is n = 
(i.e. the gate voltage V g is such that \v\ < 2). We show that in this case it is always favorable 
to slightly distort the lattice provided there is a perpendicular magnetic field 0- The distortion 
lowers the electronic energy. This energy lowering comes both from the CLL, which gives an 
essential contribution, and also from all the n < LL, which contribute in a less important way 
as we explain below. There are (2 + v)N < j > electrons in the CLL. They contribute an energy gain 




(2) 



(3) 



E, 



<n=0 — — 



N${2-\v\)M 



(4) 



1 When \v\ > 2, there is no energy gained by distortion. The distortion only occurs for v close to 
zero. 
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because when v < 0, all (2 + v)N < f ) electrons gain each an energy M and when v > 0, 2A^ 
electrons gain each an energy M but the remaining uN^ electrons loose each an energy M. This 
energy gain depends on the magnetic field through N^. In addition, the energy gain is linear in 
the out-of-plane distortion 77 because the on-site energy is proportional to the distortion, as we 
discuss below: M — Dr), where D is a proportionality constant, akin to a deformation potential 
0. The other 2(N p — N$) electrons that fill the n < LLs, also contribute to the energy lowering. 
Each of them gains a small energy compared to what an n = electron gains, as discussed in 
the preceding paragraph, but as there are many more of them, about 2(N p — N^) 2N p , we 
cannot neglect their contribution. In the tight-binding model, we find 

E n<0 = -^M\ (5) 

where the numerical factor 7 0.67. This energy gain is quadratic in the distortion, and 
therefore smaller than E n=0 at small distortion, and independent of the magnetic field. Actually, 
this term represents the full electronic energy gain for a lattice distortion under zero magnetic 
field. The distortion costs an elastic energy 

^clastic = N p G V 2 , (6) 

where the out-of-plane distortion is assumed to be small 77 <C a and G is an elastic constant, 
which - in the adiabatic approximation - is related to the ZO phonon frequency at q = by 
0Jq ~ 4G/m c , where m c is the mass of a carbon atom and ujq/2t: is the phonon frequency at the 
r point. From the measured frequency u>o/2irc ~ 800 cm -1 of the graphite out-of-plane optical 
phonon [15] . we obtain Ga 2 « m c uj 2 a 2 /A ~ 14 eV. 
The total energy is therefore 

Etat = E n=0 + E n<a + E cUstic = -N,p(2 - \u\)Dt] - -i^D 2 if + N p Gif . (7) 

The most remarkable difference in this expression with respect to the usual Peierls case is the 
presence of an energy gain linear in the lattice distortion r\ rather than proportional to rj 2 In 77 



We note that the phonon mediated (indirect) electron-electron interaction is 



el. 



Sei-ph /(h^o) ~ D 2 /G. The dimensionless electron-electron coupling constanfl A is therefore of 
order g^f'_ c \ /t. Comparing the kinetic energy gain (in the absence of a magnetic field) and the 
elastic cost in equation (0 shows that it is A — 2jD 2 / (3Gt). 

As E n< Q and -Elastic are both quadratic in the lattice distortion, we introduce a renormalized 
elastic constant G 1 = G — ^aD 2 /Hvf and write an effective elastic energy: 

^elastic + E n<0 = NpG'rj 2 . (8) 

The effect of the n < electrons is to reduce the lattice stiffness and therefore to enhance 
the distortion. We take it as an experimental fact that there is no spontaneous out-of-plane 
distortion in the absence of a perpendicular magnetic field, see also [11] . which means that 
G' > 0. This means that the Peierls instability does not occur in the absence of a magnetic 
field. For the ZO phonon, it means that - in the adiabatic approximation - its frequency 
is renormalized by the coupling to the electrons as uj 2 w uo 2 + 4(G" — G)/m c — Wq(1 — A). 
Measuring the ZO phonon mode in graphene on substrate would directly give access to G', 
and therefore provide an independent determination of the constant D through the equation 
D = y/(G - G')hv F fta. 



2 The corresponding electron-phonon coupling constant <7 c i.- P h. that would appear in a Frohlich-type 
Hamiltonian is <; i._ph. ~ Dy/h/ (m c tJo), where m c is the mass of a carbon atom. 

3 The quantity A is commonly known as the dimensionless electron-phonon coupling constant, e.g. in 
the context of the BCS theory of superconductivity, eventhough it is actually a dimensionless (phonon- 
mediated) electron-electron coupling constant [19] . 
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In order to obtain the ground state - within our mean-field approach - we minimize the 
total energy E to t, see equation (|7|), as a function of the distortion r\ and get 

1 N p 2 G" () 



from which the valley splitting follows as 



A, = ^(2-M)^«B±, (10) 

and the condensation energy is E to t = —(2 — \v\)N < f ) Dr]/2. Actually in the CLL, and only in 
that LL, the valley index is the same as the sublattice index. Therefore if at v w all electrons 
in the CLL are in one valley, they are also only on one of the two sublattices: the ground 
state can therefore be described as a commensurate A-B (or site centered) CDW. The charge 
density oscillation is quite small however: \Na — Nb\/N p w N^,/N p ~ (a/^s) 2 <§; 1, where 
I b = \Jh/{eB±_) is the magnetic length. 

The usual Peierls instability which results in a gap opening at the Fermi energy in the elec- 
tronic spectrum, is accompanied by a softening of the coupled phonon mode [19] . Approaching 
the instability from the metallic side, the phonon mode softens and goes to zero precisely at 
the transition - a phenomenon known as a giant Kohn anomaly [19] . In order to see this effect 
in the case of the magnetic field induced Peierls instability where the driving parameter is the 
magnetic field rather than the temperature, we consider broadened LL. Broadening is due to 
disorder and slightly modifies the preceding calculation for lattice distortion at weak magnetic 
field, when the valley splitting A v is smaller than the LL width Z\i mp . For example, for rect- 
angular LL - the density of states being 47V^ /A- lmp inside a LL and zero otherwise - Eq. ((4]) is 
changed into 

£ n=0 = -^M 2 = -^V, (11) 

while Eq. {2} remains unchanged since it involves totally filled LLs. The electronic energy gain 
is now proportional to rj 2 . We therefore introduce still another effective elastic constant 

G" = G'-2^-^- (12) 

in order to write the total energy as E to t = N p G"rj 2 . The renormalizcd clastic constant G" 
depends on the perpendicular magnetic field. A distortion only occurs if G" < 0, i.e. Bj_ > 
B c = 2nhG'Ai mp /3\/3ea 2 D 2 , which is always the case at large enough magnetic fields. This 
condition is precisely equivalent to requiring that the valley splitting A v = 2Drj - given by 
Eq. © with v w - be larger than the LL width A lmv . Therefore, as soon as B± is larger than 
the threshold value B c , the lattice is distorted - in other words, the ZO phonon is frozen-in - and 
the valley gap is larger than the LL width, which means that one can use the results obtained 
in the case of infinitely narrow LL. When B± < B c , the ZO phonon is strongly renormalized 
by a magnetic-field-dependent term and its frequency, in the adiabatic approximation, is 

w(B ± )*u(0)(l-2f\ l \ (13) 



where lo(0) w loq\/1 — A is the renormalized phonon frequency at zero magnetic field. The phase 
transition toward the CDW ground state occurs at B± = B Cl where oj(B±) = 0. 

Here we exhibit a microscopic mechanism providing the coupling of the electrons to the 
ZO phonon - i.e. the non-zero deformation potential D - via the interaction of the carbon 
atoms to the substrate. It is quite difficult to accurately predict the constant D and we will 
therefore only provide an order of magnitude estimate. The mechanism that we think gives 
the largest contribution results from the interaction of a single carbon atom from the graphene 



6 



Will be inserted by the editor 



sheet with the amorphous Si02 substrate treated as a dielectric continuum. The non-retarded 
Lennard-Jones interaction energy of an atom at a distance r of a dielectric wall is given by 
E-Lj(r) sa — (e r — l)(d 2 )/(e r + l)487reor 3 , where (d 2 ) is the atomic ground state expectation 
value of the squared electric dipole moment |20j . The on-site energy change resulting from the 
lattice distortion may be estimated as 

±M a E LJ (d ± rj) - £ L j(d ) a ± €r ~ ] (d2) i V (14) 

where the ± sign refers to sublattice A (+1) or B (—1), do is the average distance between the 
graphene sheet and the substrate and we assumed that 77 <C do- For a carbon atom yj (d 2 ) ~ 
4eao, where 00 is the Bohr radius, which gives Da ~ a(e r — l)e 2 a^/{e r + l)7reodQ ~ 1 to 
14 eV depending on do ~ 0.1 to 0.2 nm. Therefore, the order of magnitude of the deformation 
potential Da is 5 eV. The condition G' > that no Peierls instability occurs when B± = 
implies that Da < -J 1 Gahvp ss 9.8 eV. Therefore, in order to match the experiment [I], we 
take the plausible value Da = 7.8 eV, which gives G'a 2 w 4.2 eV. In the end, the dimensionless 
electron-phonon coupling constant A = 2^D 2 /(3Gt) « 0.8, which is quite a large value. 

The mechanism proposed here features several differences with respect to the standard 
Peierls instability [17119] : (i) the metal is two- rather than one-dimensional; (ii) the distortion 
does not change the Bravais lattice but only the point group symmetry - the inversion symmetry 
of the honeycomb crystal is broken by the instability but the triangular Bravais lattice is not 
affected by the distortion; (iii) the electron-phonon coupling does not exist in the isolated 
graphene crystal but is provided via the coupling to the substrate; (iv) the control parameter 
for the instability is the magnetic field rather than the temperature; (v) the metal (undoped 
graphene) is actually a zero gap semiconductor and has kp = and the nesting condition is 
|Q| = 0, which indeed equals 2k p. 

4 Charge density wave in graphene: direct electron-electron interactions 

Many authors have considered the effect of direct (Coulomb) electron-electron interactions in 
graphene, in particular to study the quantum Hall effect. Depending on which aspect of the 
Coulomb interaction is thought to dominate, different scenarios have been proposed. A popular 
scenario is that of quantum Hall valley and spin ferromagnetism, see e.g. |9I21I22I23I24| . This is 
the graphene version of the well-known quantum Hall ferromagnetism: the exchange interaction 
between electrons favors ferromagnetism inside a partially filled LL. Indeed, it allows one to 
reduce the interaction energy - in other words, to gain exchange energy - at no kinetic energy 
cost when the LL are assumed infinitely narrow and flat. This mechanism predicts plateaus at 
every integer filling factor in contradiction with the existing experiments. A possible way out of 
this contradiction is to consider strong disorder, which naively speaking, broadens the LL and 
therefore, by restoring the possibility of kinetic energy changes, works against ferromagnetism 
[21]. This could explain the absence of certain QH plateaus at magnetic fields up to 45 T. 

Apart from QH ferromagnetism, we are aware of only one other scenario: namely that of a 
commensurate A-B CDW ground state triggered by the magnetic field. This scenario exists in 
three different microscopic versions depending on which electron-electron interaction is believed 
to dominate. One version - the magnetic field induced Peierls instability - is the subject of the 
present paper and relies on electron-phonon coupling providing indirect interactions between 
electrons. The other two rely on Coulomb (direct) electron-electron interactions: one - the so- 
called magnetic catalysis - focusses on long-range Coulomb interactions, while the other relies 
on lattice-scale (extended Hubbard type) Coulomb interactions. All of these mechanisms lead 
to a similar ground state and therefore predict the same QH plateaus. 

4.1 Long range Coulomb interactions 

In the context of highly oriented pyrolytic graphite, Khveshchenko suggested that long range 
Coulomb interactions could also lead to a CDW ground state in the presence of a magnetic 
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field [6|. This idea stemmed from the phenomenon of dynamical chiral symmetry breaking 
originally proposed in relativistic theories of (2+l)-dimcnsional interacting Dirac fermions. 
In this scenario, the electron-electron interaction comes from the Coulomb interaction in the 
continuum limit in which only the long range part of the interaction survives. Within the CLL, 
electrons and holes pair up to form excitons, as a result of the attractive Coulomb interaction, 
and these excitonic pairs Bose condense. This leads to the creation of an excitonic gap in the 
electronic spectrum, which is similar to what we called the valley gap A v , but has a different 
magnetic-field dependence ^4 C xciton ~ s 2 /(cIb) oc \/B± 0. This mechanism only occurs in the 
presence of a magnetic field and has therefore been dubbed the "magnetic catalysis" of the 
excitonic instability [6|7j . The ground state can be described either as a Bose condensate of 
excitons inside the CLL or as an A-B CDW. 



4.2 Lattice-scale Coulomb interactions 

Several years ago, in the context of STM studies on graphite surfaces, Tchougreeff and Hoffmann 
|25l26j suggested that graphene could have a CDW ground state 0- They considered Coulomb 
interactions between electrons in an extended Hubbard model on the honeycomb lattice with 
on-site repulsion U ~ e 2 /(ea) and nearest-neighbor repulsion V ~ U, and neglected the long- 
range part of the Coulomb interaction. They showed that if the nearest-neighbor repulsion V 
was strong enough to overcome the effect of the on-site repulsion U, a CDW ground state 
should be favored. In the opposite case, a spin density wave (SDW) would result. Recently, 
several authors have shown that a magnetic field should help to trigger this CDW instability 
[9 1 8] - or the SDW depending on the precise U/V ratio [8]. In the CDW case, the corresponding- 
valley gap is linear in the magnetic field ^Hubbard ~ U(cl/Ib) ^ ^-L J us ^ as m the magnetic- 
field-dependent Peierls mechanism A v ~ (D 2 /G')(o,/Ib) 2 ■ 



5 Conclusion 

We have reviewed the proposal of a magnetic-field-induced Peierls instability in graphene lead- 
ing to a CDW ground state [5j. The driving mechanism is the electron - ZO phonon coupling, 
which is provided via the interaction of carbon atoms to the substrate. The instability cor- 
responds to the spontaneous breaking of the inversion symmetry of the graphene honeycomb 
lattice and results in Dirac electrons having an effective mass proportional to the magnetic 
field. Apart from electron-phonon interactions, there are other interactions leading to a CDW 
instability triggered by a perpendicular magnetic field |6|7|9|8j . As to the question of which 
is the dominant interaction leading to the instability in graphene, we mention two related ex- 
amples of one-dimensional metals featuring a CDW instability. On the one hand, consider the 
well-studied case of polyacetylene, which dimerizes at low temperature and has long been de- 
scribed as being the prototypical example of a Peierls instability driven by electron-phonon 
interactions. It is now known that Coulomb interactions also play a very important role in the 
instability and that electron-phonon coupling alone is not enough to quantitatively explain the 
measured effect [28]. On the other hand, in a one-dimensional metal like KCP, the dominant 
mechanism leading to the Peierls distortion has been shown to be electron-phonon interactions 
|29j . In the case of graphene on substrate, it is therefore not clear at the moment whether a 
CDW does occur and if so, which is the driving mechanism leading to such a ground state in 
the presence of a magnetic field. 



4 This excitonic gap is of the same order of magnitude as the LL spacing h\u c as the "graphene fine 
structure constant" e 2 / (AnehvF) ~ 0.5. 

Their work was actually triggered by an early experiment on a graphene sheet on the Pt(lll) 
surface |27| 
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